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QUBIT READOUT VIA CONTROLLED COHERENT TUNNELLING TO PROBE STATE 
Technical Field 

5 This invention concerns quantum computers in which the qubits are closed systems, in that 
the particle or particles are confined within the structure. In particular the invention 
concerns a closed three-site quantum particle system, in second and third aspects it 
concerns a readout method for a closed three-state quantum particle system. 

10 Background Art 

The experimental observation, manipulation and utilization of coherent quantum 
mechanical properties in solid-state systems are key technological challenges for this 
century. The importance of incoherent quantum properties has been essential for the 
15 development of microelectronics and it is hoped that coherent quantum effects will spawn 
new technologies including, but not necessarily limited to, quantum computers [34]. 

There is at present limited experimental experience in coherent solid state systems without 
optical interactions, being mainly limited to superconducting systems of qubits [35,43,33] 

20 coupled 2 qubit systems [37,27] and 2D electron gas systems [47]. Although the successes 
of the superconducting work have been impressive, there is a strong impetus to develop 
coherent technologies that are compatible with the conventional semiconductor industry, 
due to the mature manufacturing technology and scalability advantages [28,32,24,9]. 
Although a truly coherent qubit is yet to be realized in a semiconductor system, incoherent 

25 precursors to qubits have already been fabricated, for example [41,8,6] and the rapid 
progress to date suggests that truly coherent effects will be observed soon. 

An early suggested element for a quantum computer is the so-called 'charge qubit'. The 
charge qubit is where a charged particle, usually an electron or Cooper-Pair, can be in one 
30 of two distinct spatial positions. With electrical control and quantum coherence, such a 
two-state system can be used as a qubit. The concept of a charge qubit as a scalable system 
for performing quantum operations in a solid-state environment goes back to early work by 
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Ekert et al [13,2,14] and Landauer [31]. The attractiveness of such schemes is the relative 
ease of control and readout, and the obvious progression of such concepts from the 
incoherent control of conventional computing architectures. Control can be performed 
using gate electrodes, and readout via sensitive electrometers, for example single-electron 
5 transistors, SETs, or quantum point contacts. 

One important system that has been suggested for the realization of a charge qubit is the P- 
P + charge based quantum computer [24]. In this scheme, the qubit is defined by an electron 
localized to either the left or the right phosphorous ion, which constitutes two sites. This 
10 scheme has some interesting advantages over other, competing approaches, namely that 
readout is fairly easy to be achieved with single-electron transistors (SET) and gate 
operation time should be reasonably fast (~5Qps). 

A major problem with the charge scheme is the relatively high decoherence associated with 
15 distributions of charge. In fact the coupling which so readily provides the readout, is also 
responsible for the rapid dephasing. One method for ameliorating this problem has been 
suggested, namely operation in the so-called superposition basis, and preliminary 
experiments in superconducting systems are suggestive of significant improvements in 
dephasing being attainable [43]. In the superposition basis, one operates near a degeneracy 
20 point, so that the qubit is not defined by the charge being on the one site or the other, but 
rather by the symmetric or anti-symmetric combination of sites. Operation at this point has 
a greatly reduced sensitivity to noise as the potential landscape is reasonably flat. 

A further problem related to the decoherence is that population decay can occur in position 
25 space qubits on timescales which may be short compared to the measurement time needed 
for single shot readout with an electrometer, see for example Buehler et al. [7]. 

Another reason for pursuing measurements in the superposition basis derives from the need 
to perform state tomography [25] in order to characterize qubit performance. In state 
30 tomography the entire density matrix of the qubit (or more generally of the qubit system) 
after gate operation is mapped out. It requires access to non-orthogonal bases to work, and 
therefore access to both the position and superposition bases is necessary. 
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Until now, however, there has not been a natural method for performing readout of the 
superposition states and it is this problem which is the subject of this patent. For clarity, 
we will focus the following discussion of the invention to the P-P + paradigm for quantum 
computing, however it will be readily seen that the invention can be applied more 
5 generally. 

Summary of the Invention 

In a first aspect, the invention is a closed three-site quantum particle system, comprising: 
10 A first site in which the energy is controllable; 

A second site in which the energy is controllable; 

States in the first and second sites that are strongly coupled to each other by 
coherent tunnelling, and where the tunnelling rate is controllable; 

A third site in which the energy is controllable, and where the state in the third site 
15 is weakly coupled by coherent tunnelling to the first and second states, so that the third 
state is able to map out the populations of the first and second states as its energy is 
scanned with respect to the first and second states. 

A "site" could be produced by any method of confining an electron or other quantum 
20 particle, such as a dopant atom, a quantum dot, a cooper pair box, or any combination of 
these. 

The first and second states of the system may be treated as a solid-state charge qubit with 
one particle, for instance an electron, shared between the two sites. This can be operated in 
25 the superposition basis with an integrated readout using the third site as a probe state. The 
superposition basis information, which is difficult to discriminate between using an 
electrometer in a two state system, is converted using the invention into charge distribution 
information, which is relatively easy to measure. 



30 The mapping of the energy of the superposition basis constitutes a form of bias 
spectroscopy which is useful for characterizing physical parameters of a charge qubit. 
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To improve the readout mechanism a method akin to Adiabatic Fast Passage (AFP) may be 
employed. AFP is a form of charge shelving which allows close to 100% readout fidelity, 
for a review of adiabatic passage techniques, see, for example [44]. The necessary and 
sufficient requirements are a mapping of bases into different energy states which are 
5 resolvable within the linewidlh limitations of the controlled tunnelling. This allows the 
particle to remain in the third state for a length of time much greater than the measurement 
time required for the electrometer, and therefore this concept is potentially important for all 
readout schemes of qubits. 

10 The mapping requires two or more parameters to be modulated to force the system to 
remain on a specific adiabatic pathway, these parameters may be die difference between 
the probe energy and the qubit, and the tunnelling rate between the probe and the qubit. 
This method may be able to operate compatibly with re-initialisation by reversing the AFP 
trajectory so that the qubit is set into a known state. 

The energies and tunnelling rates may be controlled using gate electrodes, in particular by 
controlling the voltages on the gate electrodes. A sensitive electrometer may be provided 
for reading out the population in the third state. The electrometer may be realised by a 
single electron transistor (SET) [15] which monitors the charge in the third (weakly 
20 coupled) site, or a quantum point contact. 

Aside from the possibility of increased decoherence times, there are other advantages to 
superposition state readout. State tomography is necessary in order to properly characterize 
qubit operation. Tomography requires measurements in more than one basis, and therefore 
25 access to alternate bases is essential [25]. The invention is useful for combined 
superposition and position basis measurements. 

These concepts are applicable to readout from a range of quantum computers. The 
quantum computer may use charge qubits, such as semiconductor donors, semiconductor 
30 quantum dots or superconducting Cooper pairs. 



Alternatively, the system may involve a solid-state spin qubit on the first site, with a 
reference spin on the second site, and the third site may be used as a probe site. The 
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relative spins of the first and second sites (qubit and reference) are impossible to 
discriminate using an electrometer, but they can be converted using the invention into 
charge distribution information, which is relatively easy to measure. 

5 The mapping of the energy of the spin states constitutes a form of bias spectroscopy, which 
may be useful for characterizing physical parameters of the spin qubit and the two spin 
system. 

To improve the readout mechanism a method akin to Adiabatic Fast Passage (AFP) may be 
10 employed. AFP is a form of charge shelving which allows close to 100% readout fidelity, 
for a review of adiabatic passage techniques, see, for example [44]. The necessary and 
sufficient requirements are a mapping of bases into different energy states which are 
resolvable within the linewidth limitations of the controlled tunnelling. The mapping 
requires two or more parameters to be modulated to force the system to remain on a 
15 specific adiabatic pathway, these parameters may be the difference between the probe 
energy and the qubit, and the tunnelling rate between the probe and the reference states. 
This method may be able to operate compatibly with re-initialisation by reversing the AFP 
trajectory so that the qubit is set into a known state. 

20 The energies and tunnelling rates may be controlled using gate electrodes, in particular by 
controlling the voltages on the gate electrodes. A sensitive electrometer may be provided 
for reading out the population in the third state. The electrometer may be realised by a 
single electron transistor (SET) [15] which monitors the charge in the third (weakly 
coupled) site, or a quantum point contact. 

25 

These concepts are applicable to readout from a range of quantum computers. The quantum 
computer may use spin qubits, such as nuclear spin or electron spin. 

The invention may also be useful for other physical realizations of qubits, such as those 
30 involving flux quanta or 2D electron gas, in which a third state is introduced for readout. 

In a second aspect, the invention is a readout method for a closed three-state quantum 
particle system, comprising the following steps: 
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Controlling the energy of a first site; 
Controlling the energy of a second site; 

Controlling the tunnelling rate between the first and second sites, which are strongly 
coupled to each other by coherent tunnelling; 
5 Controlling the energy of a third site, the state in the third site being weakly coupled 

by coherent tunnelling to the first and second states, so the third state is able to map out the 
populations of the first and second states as its energy is scanned with respect to the first 
and second states. 

10 To improve the readout mechanism a method akin to Adiabatic Fast Passage (AFP) may be 
employed. AFP is a form of charge shelving which allows close to 100% readout fidelity, 
for a review of adiabatic passage techniques, see, for example [44]. The necessary and 
sufficient requirements are a mapping of bases into different energy states which are 
resolvable within the linewidth limitations of the controlled tunnelling. This allows the 

1 5 particle to remain in the third state for a length of time much greater than the measurement 
time required for the electrometer, and therefore this concept is potentially important for all 
readout schemes of qubits. 

The mapping requires two or more parameters to be modulated to force the system to 
20 remain on a specific adiabatic pathway, these parameters may be the difference between 
the probe energy and the qubit, and the tunnelling rate between the probe and the qubit or 
reference state. This method may be able to operate compatibly with re-initialisation by 
reversing the AFP trajectory so that the qubit is set into a known state. 

25 The energies and tunnelling rates may be controlled using gate electrodes, in particular by 
controlling the voltages on the gate electrodes. A sensitive electrometer may be provided 
for reading out the population in the third state. The electrometer may be realised by a 
single electron transistor (SET) [15] which monitors the charge in the third (weakly 
coupled) site, or a quantum point contact. 

30 

In a third aspect the invention is a readout method for a closed three-state quantum particle 
system, comprising the following steps: 
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Mapping bases into different energy states which are resolvable within the 
linewidth limitations of controlled tunnelling to a probe state, where the mapping requires 
two or more parameters to be modulated to force the system to remain on a specific 
adiabatic pathway. The parameters may be the difference between the probe energy and the 
5 qubit, and the tunnelling rate between the probe and the qubit or reference state. 

This method may be able to operate compatibly with re-initialisation by reversing the AFP 
trajectory so that the qubit is set into a known state. 

10 The energies and tunnelling rates may be controlled using gate electrodes, in particular by 
controlling the voltages on the gate electrodes. A sensitive electrometer may be provided 
for reading out the population in the third state. The electrometer may be realised by a 
single electron transistor (SET) [15] which monitors the charge in the third (weakly 
coupled) site, or a quantum point contact. 

15 

Brief Description of the Drawings 

Examples of the invention will now be described with reference to the accompanying 
drawings, in which: 

20 

Fig. 1 (a) is a schematic diagram of a qubit, probe and readout gate; Figure 1(b) is a triple 
well energy diagram for the qubit and probe. 

Fig. 2 is a state energy diagram for the qubit of Fig. 1 . 

25 

Fig. 3 (a) is a graph showing the variation of p u (black solid line), pn- (black broken line) 
and p pp (grey solid line) as a function of time (in units of n/Q/ r ) for A pi /Q fr = 2 . Fig. 3 
0>Mg) are density plots showing the populations of p u in (b) and (e), p^ in (c) and (f), 
and p^ in (d) and (g), as a function of time and A pt (in units of Q /r ) for £/=£^0, 
30 £1^ =Q^. /20 and r = £\ /100. 

Fig. 4 (a)-(c) are density plots showing the populations with /<0)=(1/2)[1;-1;0;-1;1;0;0;0;0] 
(initialized in the anti-symmetric state), (d)-(f) are density plots showing the populations 
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for p(0)=(l/2)[l ; 1 ;0; 1 ; 1 ;0;0;0;0] (initialized in the symmetric state). The density plots show 
the populations p u in (b) and (e), PtT in (c) and (f), and Ppp in (d) and (g), as a function of 
time and A pl (in units of Q fr ) for £>=£^=0, Cl^ = Cl tr 120 and T = Q /r /100 . 

5 Fig. 5 (a)-(c) are density plots of populations as a function of time and A pl with p„ (0) = 1 
(initialized in |/». p u {t) in (a), Prr {t) in (b) and p pp (t) in (c) and (d), a smaller colour 
axis is used in (d). 

Fig. 6 (a)-(c) are density plots of populations as a function of time and A pl with p„(0) = 1 
10 (initialized in |/». p u {t) in (a), Prr {t) in (b) and p pp {t) in (c). The difference between 
these traces and those in Fig. 5 (a) to (c) is the enhanced Cl^ = 0.25n fr . 

Fig. 7 (a) is a graph of as a function of time for the AFP style population transfer. Fig. 

7 (b) is a graph of the energy eigenstates as a function of time for the AFP style population 
15 transfer. 

Fig. 8 (a) is a graph of populations as a function of time for the AFP trajectory in Eq. 3 for 
the anti-symmetric state, (b) is a graph for the symmetric state, (c) is a graph of a mixed 
state with equal initial populations in |/> and |r>, (d) is a graph of the results for the anti- 
20 symmetric state/>„(0) = l (|/» and (e) is a graph of the results for the symmetric state 
Ar(°) = 1 In each case the solid black line correspond to p u {t), the dashed black line 
to p„ (t), and the solid gray line to p pp (t) . 

Fig. 9 is a graph showing the readout fidelity (solid black) and population transfer 
25 probabilities from the symmetric (solid gray) and anti-symmetric (dashed black) states as a 
function of for fi™ = 0.8 

Fig. 10 is a schematic diagram for a proposed system for tomographic readout of the charge 
qubit. 

30 

Fig. 11 (a) is a schematic diagram of a simple scheme for spectroscopic readout. Fig. 11 
(b) is a schematic diagram of a more complex scheme where an extra gate is used to tune 
the coupling. Fig. 11 (c) is a schematic diagram for a more complex scheme where 
coupling is controlled via flux threaded through a Cooper-pair box. 
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Fig. 12 (a) is a schematic diagram of an arrangement for superposition state readout which 
will be applicable to the readout of spin qubits. Fig. 12 (b) is a triple well energy diagram. 

Fig. 13 (a) is a graph of the eigenvalues for a two-spin, three-well case as a function of the 
5 energy of the probe state, Ep. Fig. 13 (b), (c) and (d) are transient bias spectroscopy 
readouts showing readout S as a function of time in units of nfi/h as a function of Ep/Ji,, 
with (b) showing any of the triplet states, (c) showing the singlet state, and (d) showing a 
superposition of the singlet and one of the triplet states. 

10 Fig. 14 is a series of graphs against time. Fig. 14(a) shows the AFP pulse. The left axis is 
Qq/Jir, corresponding to the Gaussian variation in the rp tunnelling rate, the right axis is 
Ep/Jt corresponding to the linear ramp from high to low of the probe. Fig. 14(b) is a time 
varying eigenspectrum depicting evolution through the AFP sweep. Note that the form of 
the eigenspectrum is slightly modified from those shown earlier as £lrp is now varying as a 

15 function of time. Fig. 14(c) is a readout statistic for the triplet states (uppermost traces), 
superposition of singlet and triplet states (central traces) and singlet state. 

Best Modes of the Invention 

20 Referring first to Fig. 1(a), a three-donor system comprises three ionized phosphorus 
donors 10, 11 and 12 sharing a single electron (not shown). The strongly coupled qubit is 
defined by the donors 10 and 11 labelled |/> and |r>. The weakly coupled probe 12 labelled 
|p> provides the readout. We follow the gate notation used in Refs. [24]. 

25 The energies of each single-electron state are controlled using shift gates 15, 16 and 17, 
labelled Si, S r , S p , and the energies of these states are £/, E r , Ep. 

The coherent tunnelling rates, fi> and Cl^ are controlled by barrier gates 20 and 21 labelled 
Bif and Bjj, respectively, which control the barrier height (and hence the coherent 
30 tunnelling rates). There is assumed to be no tunnelling between 10 and 12. The B 
dependence is highlighted as a dynamically tunable variable, but the static properties 
dependent on fabrication processes have been dropped. This follows the spirit of Gurvitz's 
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[19] treatment for the two coupled-dot system and accords with the work of Renzoni and 
Brandes [38] in a triple-well system. 

An SET 25 reads out the electron on the probe donor 12. 

5 

The natural basis for the qubit will be the superposition basis, we write: 
the symmetric state as |S) = (l/V2)(jr)+|/)) and 

the anti-symmetric state as \AS) = (l/V2)(j r) -| /)). 

10 To readout, the energy on the probe 12, E p , is varied using gate 17 and the population in 12 
is monitored with the SET 25. There will be nonzero populations in 12 only when E p is 
degenerate with either of the superposition states E\ or E r because these do not correspond 
to eigenstates of the strongly-coupled system. Similar bias spectroscopy is seen in 2DEG 
systems [47]. In optical Autler-Townes (AT) [1] experiments, it is usual to monitor the 

15 response of a weak probe field, which is proportional to the coherence pr P , but in this case 
we are required to monitor the population . In Fig. 1(b) the triple well diagram shows a 
possible configuration where Ej=E r so that symmetric |5> and antisymmetric \AS) states are 
equally separated from Ej=E r . Ep has been tuned to E AS so that tunnelling from \AS) to \p) 
occurs. 

20 

In order to proceed, we will investigate the Hamiltonian to show the robustness of the 

superposition state to charge noise, and then numerically solve the density matrix equations 

of motion to show the transient and spectroscopic details. This is similar to the approach in 

Ref. [16]. We write down the interaction Hamiltonian (all energies expressed relative to 
25 Ei): 

// = Ajr)<r| + A>>^ (1) 
where A^ = E a —E l9 a = r,p. This Hamiltonian is identical to the well known three-level 
atom Hamiltonian in the rotating wave approximation from optics [40]. 
30 Eq. 1 may be more conveniently expressed in matrix form as: 
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0 0 

o an a„ 



Fig. 2 shows the energy levels of the |/>-|r> qubit system as a function of £/-£ r , which is a 
linear function ofS r Si for a given Q, r (with all quantities normalized to Q fr ). Note that 

5 the eigenstates are presented unnormalized in the figure. The solid lines are the energy 
eigenstates, the dashed lines correspond to the non-interacting (Q /r =0) states. As the 

degeneracy point is approached, the eigenstates adiabatically evolve from the position 
states to the superposition states when coherent tunnelling is allowed by the barrier gates 
20 and 21. The lower lying state at the degeneracy point is the symmetric state, and the 
10 higher, the anti-symmetric state. 

"Adiabatic" is a word describing how evolution occurs. Adiabatic evolution is where a 
system is allowed to evolve so that it stays in the 'steady state* at all times - the system is 
always in equilibrium with the forces driving it. Usually, one applies slowly varying stimuli 
15 to achieve adiabatic evolution. Non-adiabatic (or sometimes diabatic) is the opposite of 
this, where a system is excited by a stimulus that can cause an abrupt change of state, and 
then the system responds causing c ringing\ 

The experiment directly analogous to the Nakamura experiment [35] in a charge qubit 
20 would be to initialize the system in state |/> (far left, start point) by setting We next 

shift the system nonadiabatically to the degeneracy point, £/=£>, allow coherent evolution, 
and then nonadiabatically return to the start point (E F »E0. As a function of the time spent 

at the degeneracy point, the final population in |r> will show characteristic oscillations, the 
signature of coherent evolution. 

25 

Noise can be viewed as a perturbation of the shift gates 15, 16 and 17. It is immediately 
clear from Fig. 2 that the system will be far more sensitive to gate noise well away from the 
degeneracy points. The rate of change of the energy level with gate voltage is minimized at 
the degeneracy point. We therefore would expect the system to be far more robust against 
30 decoherence mechanisms at the degeneracy point, than when well away from degeneracy, 
again, suggestive of the importance of superposition state operation and readout. 
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We now turn to the problem of solving the transient spectroscopy by investigating the 
density matrix equations of motion 

5 p = -L[H 9 p]+L 

where p is the density matrix, and Lis an operator describing dephasing. 

The T2 processes are modelled by introducing a dephasing rate T which is assumed to 
couple to all of the coherences, the T\ processes are described by rates of incoherent 
10 population transfer [42] 

r = y A <* 



lafi l-expiA^/kT) 

where k is Boltzmann's constant, T is the temperature and x a p is some ra te related to the 
tunnelling probability between states a and p. is the rate of population transfer from (3 

15 to cl Note that we have assumed no direct population transfer between states |/> and \p) by 
any means. 



From here we may write down the density matrix equations of motion, which are: 
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Pi, = &lr (Pr,-Pu)~ r rt p u +T lr p rr , 
Pu=^-P t r +n i r(Prr-P«)-nr P P l p]-r Plr , 

Pip = Pip + V-lrPrp ~ torpP^ ~ rpl P > 

Prr=i[^lr(P l r-Pr,) + ^r P (P P r ~Pj]+K,Pu ^pr)prr + ^r P Ppp> 

Prp=\^-Prp +n rp (p pp - Prr )+il lr p lp ^-r Prp , 

Ppp = ' n r P (Prp -P^^prPrr- T rpP Pp , 

* 

Pofl ~ Pafi > 

1 = Pll+ Prr + Ppp- 

(2) 

We solve Eqs. 2 by converting to a set of matrix equations. 

5 

X=(P+L)Y 

where 



Pit 






Plr 




Pa 


P'P 




Plr 


Prl 




Pip 

Prl 


Prr 


,Y = 


Prr 


Prp 




Prp 


Ppl 




Ppl 


Ppr 




Ppr 


Ppp 






1 







10 and 
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Pn 
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n fr o 
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0 








-Q fr 0 
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rp > 
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-£} A /ft 


-A,,/* 




0 








0 




0 



and i> 2 = P 2l = ifi fr I(3) , P u =P 32 =/Q r/) I(3), where 7(3) is the 3x3 identity matrix, 
^13=^31=0 and P 4>1 =[1,0,0], P 4>2 =[0,1,0], P 43 = [0,0,1] . We also have 



£ = - 
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10 We study these equations using two approaches. Firstly we employ a spectroscopic 
approach to study p pp as a function of time and detuning, A p/ . Secondly we map A / and 
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along a specific trajectory to map population from one of the superposition states to 
state |p) using an Adiabatic Fast Passage (AFP)-like method. 

Results and Discussion 

5 

In this section we present results obtained from the investigations of the density matrix 

equations of motion, Eqs. 2 looking at the transient response. These results will be shown 

to be analogous to similar optical cases, such as those presented in [16]. The two following 
sections deal with the time-dependent spectroscopy, i.e. for constant A p/ , and the use of the 

10 third dot as a readout for a superposition based quantum computer, and therefore shows the 
result of performing AFP oh the dot. 

Time-dependent spectroscopy 

15 Almost all systems that undergo non-adiabatic evolution exhibit some kind of ringing, and 
the quantum computer systems are no exception. The periodic oscillations presented in 
Fig. 3 are population oscillations due to the system's having been artificially initialised in a 
state that is not an eigenstate of the system (and therefore not in the steady state). 
We have numerically integrated Eqs. 2 for various cases with p u (0) = 1 so as to highlight 

20 all of the spectroscopic details. We also set Ef=E^=0 9 Cl^ = Q /r /20, T = Q, r / 1 00 and 
Zap = 0 . In Fig. 3: (a)-(d) show populations for p a (0) = 1 , and (e)-(g) show populations 
for p rr (0) = 1 . Notice the dominant oscillatory behavior in p u (t) and p^if) and the AT 
doublet-like feature in p pp . Fig. 3 (a) shows the time dependent populations when E p =Eas- 

The dominant feature is the coherent population oscillations between / and r. There is also 
25 a steady buildup of population in p which is our measurement signal. In Figs. 3 (b)-(d) we 
show piu p rr , and p pp respectively, as a function of time and A^. Again, the dominant 

behavior is the coherent oscillation between / and r, however when = ±Q /r , resonant 

tunnelling into p can occur, leading to the appearance of a doublet in p pp similar to the 

optical AT doublet. This population signal shows the bias spectroscopy allowing access to 

30 the superposition state dynamics. Analogous to Figs. 3 (b)-(d), in Figs. 3 (e)-(g) we show 

the dynamics when the system is initialized in state |r>. The fact that both sets of dynamics 
show very similar results for p pp can be easily understood by realising that our 
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measurement basis is the superposition basis, and \l) and |r> can each be expressed as equal 
superpositions of symmetric and anti-symmetric states. 

In order to prove that we are indeed able to discriminate between the symmetric and anti- 
5 symmetric basis, we begin our simulations in either of these two states and observe the 
results. To initialize in the symmetric state, we choose p(0)=(l/2)[l;l;0;l;l;0;0;0;0], and 
to initialize in the anti-symmetric state we choose, p(0)=(l/2)[l;-l;0;-l;l;0;0;0;0]. The 
obtained pseudo-colour plots are presented in Figs. 4 (a)-(c) (anti-symmetric state) and 4 
(d)-(f) (symmetric state). The desired signals are clearly seen by the nonzero p pp at the 

10 energies corresponding to the energies of the symmetric and anti-symmetric states. This 
shows that this kind of triple-dot spectroscopy is (at least in the ideal case) able to 
discriminate between the superposition bases. The presence of a small signal at the energy 
corresponding to the undesired superposition state is due to dephasing effects which take 
population from one superposition state to the other. 

15 

The oscillatory behavior characteristic of Fig. 3 is no longer present as the initial states are 
eigenstates of the strongly coupled Hamiltonian, Large population signals in p pp are only 

observed when resonant with the appropriate state. 

20 Including nonzero T\ processes complicates matters as is evinced by Fig. 5 which was 
generated using the above parameters, Ef=E^=0, and kT/(fl£l lr ) = 0.002, x lr =Zrp = 0.1. 
Figs. 5 (a), (b) and (c) shows p„, Prr and p pp respectively. Fig. 5 (d) shows p pp on a 

smaller colour axis to highlight the population transfer from the anti-symmetric state. Note 
that the dynamics for large negative A pI are dominated by the incoherent relaxation 

25 processes pumping population into state |p>. Such dynamics are clearly unsuitable for a 
quantum coherent readout. Just considering the space for A pl > 0 , however, shows that the 
behavior is qualitatively very similar to that found without T\ processes. 

The T\ dephasing processes have little effect on the resonant |/>-|r> system (qubit) except 
30 when A pl is at modest negative values, indicating that the system will relax (incoherently) 

to this state. The fine scale used in (d), however, shows that a spectroscopic readout will 
still be possible in the presence of T\ processes (albeit with reduced signal), provided we 
maintain A ,)0. 
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The above plots in Figs. 3 to 5 show the effect of a weak probe. However if the coupling to 
state \p) becomes significant, then the simple dynamics of a strongly coupled qubit probed 
via a weakly coupled dot become insufficient and true three-state dynamics need to be 
considered. To illustrate these three-state dynamics, we present in Fig. 6 a set of plots 
5 analogous to Fig. 5, but where H^/Q,, has been increased from 0.05 to 0.25. The 
combined effects of T\ and strong coupling have yielded a very complex system, which 
would appear to be undesirable for the operation of a quantum computer. 

Readout via Adiabatic Fast Passage (AFP) 

10 

Usually, one applies slowly varying stimuli to achieve adiabatic evolution, and indeed in 
most of the literature, adiabatic is taken to mean this; however it is not necessary. All that 
is required for adiabaticity is to maintain the steady state. In order to achieve adiabatic 
evolution on fast time-scales, one can employ adiabatic fast passage techniques, see 
15 Vitanov et al. [44] for a thorough review of these techniques in the optical/RF regimes. In 
these techniques a specific trajectory through phase space is taken by the system, according 
to controlled stimuli that force the system to stay on the specific adiabatic pathway. 

In general, two or more parameters need to be modulated together to achieve adiabatic fast 
20 passage, in contrast to simple adiabatic evolution where only one external parameter needs 
to be varied. The advantages of adiabatic fast passage techniques include usually being 
highly robust to errors in the applied stimuli, compared with nonadiabatic pulses. They are 
relatively fast, typically requiring of order 10 coherent oscillation periods for transfer of 
population (compared with nonadiabatic transfer which is typically of order 1 oscillation 
25 period but susceptible to noise, and adiabatic evolution, which may require of order 100 
oscillation periods). Adiabatic fast passage pulses require pulse shapes that are easily 
realized experimentally: typically exponential, or in our case cosinusoidal pulses which are 
attainable with most standard generators, contrasted with the square pulses required for 
non-adiabatic evolution which would need unfeasibly fast rise times in our experiments. 
30 The trade-off is a slight increase in complexity because of the need to modulate more than 
one parameter simultaneously, although this should be relatively easily attained and indeed 
are required for full quantum computation in our system in any case. 
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Rather than restrict ourselves to setting the detuning, A pt and reading the population which 
tunnels into state |p> as is done above, we can drive the population into state \p) (charge- 
shelving) from one of the superposition states by tunnelling control similar to AFP. The 
advantages of this are that p pp can be adiabatically driven to a large value (approaching 

5 unity) in a short time (typically a few £\~ l ), with high robustness against gate errors. 
Because of the energy dependent T\ decay processes, it is most useful to perform AFP 
between state \p) and the most energetic of the two superposition states, in this case the 
anti-symmetric combination. One needs to take some care in choosing the trajectory along 
which to modulate state \p) in this case, due to the fact that there will be off-resonant 
10 interactions. Our results are promising and suggest that optimization will be possible. 

The energy of |p> is swept through resonance with the anti-symmetric superposition state, 
with the tunnelling rate pulsed. The trajectory taken by state |p> is governed by 
controlling both A p/ and Q. rp and for the traces presented in Figs. 7 they were: 

15 

A /p =2Q /r (l-;/0 

cy™** (3) 
^=-f-[l-cos(2^/^)l 

where = 0A€l u , = 0 and = 5a/a fr . In order to make this more explicit, Fig. 
7 (a) shows Cl rp (t), note that in keeping with conventional AFP schemes, the scheme is 
fairly insensitive to the exact form of Q^; and (b) is a diagram showing the energy levels 

as a function of time for the trajectory used. Note that we have chosen a cosine 
dependence for Cl rp rather than the more usual exponential dependence to ensure that 
Cl rp = 0 at the beginning and end of the pulse. 

The maximum tunnelling rate is when E p =E^s ensuring adiabatic population transfer. 

Fig. 8 (a) shows the populations for the qubit being initially prepared in \AS} with initial 
conditions /<0)=(1/2)[1;-1;0;-1;1;0;0;0;0]. After AFP most of the population has been 
driven into |p>. Similarly Fig. 8 (b) shows the effect of the AFP trajectory on \S) with initial 
conditions /<0)=(1/2)[1;1;0;1;1;0;0;0;0;|. Note that there is minimal population transfer. It 
30 is important to realize that the superposition states contain information about the coherence 
in the system, and so it is to be expected that there will be a different signature of 



20 



25 
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incoherent states. In Fig. 8 (c) we show the results for a system initialized in a mixed state 
of |/> and |r) with initial conditions p(0)Kl/2)[l;0;0;0;l;0;0;0;0], and in Figs. 8 (d) and (e) 
we show results for initialization in states [/> and |r> respectively. In each case the solid 
black line correspond to p a (f) 9 the dashed black line to p^t), and the solid gray line to 
5 Ppp(?) • In all cases, the measurement result at the end of the AFP trajectory is p pp « 0.5 . 
This result is hardly surprising, corresponding as it does to performing a projective 
measurement onto the superposition basis, and having equal probability of the system being 
in either the symmetric or anti-symmetric state. Note that the different initial conditions 
between the mixed states and states \l) and |r) suggest a tomographic scheme. 

10 

The readout fidelity of this method can be measured by comparing the final values of p pp 
for the different basis states, i.e. 

p P PP (^niax)~ P pp max) 

15 We have calculated the readout fidelity as a function of for other parameters constant 
as above (although T\ processes turned off) which is presented in Fig. 9. As expected, the 
readout fidelity is maximized when » €l rp . This is because off resonant interactions 

between the symmetric state and |p> are minimized. In the case considered, the maximum 
fidelity attained was around 90%. Although this number is far short of that required for 
20 scalable quantum computing, it is certainly possible to optimize this fidelity, which should 
enable arbitrarily high fidelities to be obtained. 

Notice in Fig. 9 that there is a small region where the fidelity is negative, this is because, 
there is in fact stronger coupling to the undesired state than the desired state. 

25 

Note also that an important aspect of measurement of a qubit is the ability to re-initialize 
the qubit after measurement. The AFP scheme presented here is entirely reversible, and 
hence following readout it is straightforward to pump the electron back from the probe 
state into the anti-symmetric state with equally high fidelity as that obtained when 
30 performing the readout. 



Other Systems 
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In this section we present some examples of methods of implementing similar schemes to 
the readout mechanism described above, more applicable to systems other than the P-P + 
charge-qubit. We first discuss tomographic style readout, followed by mechanisms for 
observing similar effects in superconducting experiments, and then applicability to spin- 
5 based qubits. As mentioned earlier, however, these are examples of what is essentially a 
highly general scheme, and experts in the field will be able to perceive the generality and 
applicability of the concepts. 

Tomography 

10 

Tomographic measurements, essential for proper characterization of quantum 
computational performance, require measurements in more than one basis [25]. The access 
to the superposition basis afforded by our measurement scheme is therefore extremely 
advantageous to the goal of tomography of the charge qubit. In Fig. 10 we show one 
15 scheme which is suggestive of a tomographic method, the arrangement is similar to Fig. 1 
with the addition of a second SET 35 on the left side. The position basis can be monitored 
with the leftmost SET 35, and the superposition basis monitored with the rightmost SET 25 
using AFP. Performing multiple measurements in alternate bases will enable full 
tomographic readout and density matrix reconstruction. 

20 

Reconstruction of the density matrix will require many measurements in the alternate 
bases. In Table 1 we show the results seen on each SET for an idealized experiment 
averaged over many experiments (where |MS> is the mixed state, defined by the density 
matrix p=(l/2)[l;0;0;0;l ;0;0;0;0]). The measurements presented are with the AFP scheme, 
25 measurements in the position basis (AFP |f>(/|) and superposition basis (AFP |p)<p|), and 
without the AFP (No AFP |/>(/|). The fact that each of the five measured states give 
different measurement signatures demonstrates the feasibility of the tomographic scheme. 



Initial State 


AFP|/></i 


AFPbXd 


NoAFPMfl 


\l> 


0.25 


0.5 


1 


\r) 


0.25 


0.5 


0 


15) 


0.5 


0 


0.5 


\AS) 


0 


1 


0.5 


|MS> 


0.25 


0.5 


0.25 
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Table 1: Table showing measurement outcomes for the double SET system in Fig. 10. 
Measurements are.to be performed both with and without AFP. If AFP performed, then 
both SETs are read out, otherwise only the leftmost SET is monitored. 

5 

Superconducting systems 

Superconducting systems are an obvious choice for systems which exhibit coherent 
quantum effects. There are several different flavours of superconducting qubits, and we 

10 will not discuss all of them, rather concentrating on one example. The generality of these 
concepts, however, should be obvious to one experienced in the art. The Cooper-pair box 
system is a particularly interesting system, and historically has been the earliest solid state 
system to exhibit the potential for scalable qubit operation (see for example Ref. [35]). The 
schemes presented in this work are certainly applicable to the Cooper-pair box system, and 

15 this system may therefore also benefit from the increased lifetimes associated with the 
superposition state. 

In Fig. 11 we present three suggestions for superconducting circuits which would be 
appropriate for the spectroscopic and AFP style measurements presented. These 
20 superconducting circuits are based on Cooper-pair boxes. The open circles 50 correspond 
to low resistivity tunnel junctions, the closed circle 60 to a high resistivity tunnel junction. 
70 labelled B indicates a magnetic field threaded through the loop perpendicular to the 
page. 

25 Fig. 1 1 (a) is a simple scheme which will only be capable of spectroscopic readout as the 
high resistivity tunnel junction 60 will not provide a tunable coupling. 

Fig. 1 1 (b) and (c) illustrate alternative methods of performing the required AFP coupling 
schemes by tuning the coupling to the final area. Gate 80 labelled G in (b) is an extra gate 
30 used to tune the coupling whereas in (c) the coupling £1^ is controlled via the flux 70 

threaded through the Cooper-pair box. 

Note that these suggestions are by no means exhaustive, they are merely presented by way 
of example. 

35 
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Spin readout 

Because of the interest in spin based qubits, there have been many proposals for methods to 
readout the state of a single spin. One example is the charge to spin conversion proposed 

5 by Kane [28]. That scheme (proposed for a phosphorus in silicon qubit) involves 
application of external gate biases to allow spin dependent tunnelling of an electron from 
one donor onto a second donor, creating a two electron state, the D" state. This charge 
difference could then be measured using an SET. However, it is unclear if it will be 
possible to stabilize the D" state for long enough to permit a successful readout. Another 

10 scheme for readout of electron spins comes from Loss and DiVincenzo [32] who propose 
two methods for spin readout in a quantum dot system, one by switchable tunnelling into a 
cooled paramagnetic dot, the magnetization of which is sensitive to the spin of the electron 
tunnelled, the other being the use of a 'spin valve 5 (see for example Konig and Martinek 
[30]. We are not aware of any single-shot spin measurements using a spin-valve, nor are we 

15 aware of experiments towards realizing the paramagnetic dot. Magnetic resonance force 
microscopy (MRFM) has been suggested as a readout mechanism, and the obtained 
sensitivities may well reach the single-spin level in the near future [10]. However we do 
not believe that such a method will constitute an effective and scalable readout mechanism 
for a spin based quantum computer. Optical readout has also been proposed (for example 

20 the resonant micro-Raman scheme of Koiler et al [29]). Although certainly useful for bulk 
characterization, optical spot size limitations and optical collection efficiency will hamper 
effectiveness for single-shot quantum computing applications. 

We propose a modification of the superposition state readout which will be applicable to 
25 the readout of spin qubits similar to that depicted in Fig. 12 where there are three 
phosphorus donors 10, 1 1 and 12, the rightmost 12 being ionized and acting as the |p> state. 
In the figure we have shown the case in the absence of a magnetic field, so that the triplet 
states are degenerate, however with a magnetic field applied, this degeneracy will be lifted. 
If the separation of the states exceeds the effective transfer linewidth, then all of the 
30 different states will be resolvable using bias spectroscopic methods, however in order to 
use an AFP approach, it will be necessary to fine tune the parameters in order to achieve 
high fidelity readout. The spins of the phosphorus donors 10 and 11 (election or nuclear 
spins) will couple via the exchange interaction and there will be an energy difference 
between the singlet and the triplet states. In particular the two spin system on donors / and r 
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will couple via the exchange interaction to realize a singlet \S) = ^^((t>l)- |4T)) or 

triplet |7) state. The exchange interaction lowers the energy of the singlet state relative to 
the triple, and so energy dependant tunnelling (bias spectroscopy) will be possible as 
before, with tunnelling from r to p occurring when Ep=Ef or E p =E$ as appropriate, 

5 yielding a single shot spin readout, compatible with conventional RF-SETs. This all- 
electrical readout mechanism is therefore promising and highly applicable to the single- 
shot readout required for a quantum computer and is potentially superior to all other 
proposed spin readout mechanisms for this purpose. 

10 We introduce here the full Hamiltonian in second quantization form for the two electron, 
three site problem (the form is different from the form used in Eq. (1). 

H = 5[Z*A!A« +*&MaKa +blbj+m„(b: ta b p , a +^.aJ 

+ 42/,(&iWW -blb ia b)-b fa ^\blb^b) a b f ^ +2^X^*fe^t -bib*) 



(4) 



15 where i and j are summed over all sites; a is summed over f ,>l with a =4<,T ;b ia is the 
annihilation operator for an electron on site i with spin a; B is the magnet field; E x is the 
electrostatic energy of an electron on site i; Jy is the exchange energy between sites ij with 
Jn=0; and we define B*=2g* mbB as the Zeeman energy splitting. 

20 The first point of attack for studying a Hamiltonian like Eq. (4) is to perform an eigenvalue 
analysis to identify the states and their migration through phase space. Anti-crossings 
between levels indicate interactions, and crossings indicate non-interacting levels. 

Fig. 13 (a) is a graph of the eigenvalues for a two-spin, three- well case as a function of the 
25 energy of the probe state, Ep, for constant Zeeman splitting and exchange interaction. 
Vertical lines correspond to eigenstates without any component of |p> in them. The 
diagonal states have some |p) and are split into two submanifolds, corresponding to 
symmetric (left) and anti-symmetric (right) combinations of |1> and |r). 
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Transients are obtained by solving the density matrix equations of motion. For simplicity 
we have not included dephasing processes, although these rates are expected to be 
compatible with those assumed in charge and spin based qubit schemes (note that we 
require only Cl^ to be greater than the charge coherence time, and the system is only 
5 vulnerable to charge relaxation during this procedure, this requirement is unchanged from 
the initial proposal). 

The density matrix equations of motion are 
10 P = ~[H,p] 

which we solve numerically for various initial conditions. A realistic experiment to show 
these results would comprise (i) initialisation in a known state, (ii) evolution for length r , 
(iii) measurement. The time axis in the following discussion is therefore indicative of the 
15 evolution as a function of the Tiold 1 time, r . 

Fig. 13 (b), (c) and (d) are transient bias spectroscopy readouts showing readout S as a a 
function of FJi^ and time (in units of nh/h for B*/Jir=0.4, Clrpfh = 0.1, Ei = E,- = 0 and 
various initial conditions, with (b) showing any of the triplet states, (c) showing the singlet 
20 state, and (d) showing a superposition of the singlet and one of the triplet states. 

As with the need to operate with Cl tr > J lr (which is physically realistic) observed 
previously in the eigenvalue and transient analysis, so to with Adiabatic Fast Passage 
(AFP) we need to shift the overall working point along the E p axis, but there are no further 
25 modifications. The robustness and tolerance of the scheme is illustrated in Fig. 14, which 
shows AFP pulses, time varying eigenspectra and readout fidelity for Q /r = J u . As Qi/Jir 

increases above unity, the fidelity of the transfer is found to improve. 

The form of the AFP sweep has been discussed earlier, this time we choose: 

30 
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10 



15 



where = 0.257,^ = 0,<r =/ max /8 the standard deviation of the pulse, and 

'max = 60^/ •/ . Note the only qualitative difference between this form and that presented 

5 earlier is the shift in the absolute position of E p (the change from cosinusoidal to 
exponential form is trivial). 



Discussing first Fig. 14, the top trace (a) shows the pulsed rp coupling, IJ U (left axis) 
and probe bias Ep/J^ (right axis) as functions of time. The corresponding time varying 
eigenspectrum is presented in the middle trace (b). The presence of the extra peaks at the 
top of the eigenspectrum is due to the other submanifold of states, however as can be seen 
they do not appear to interfere with the passage. The bottom trace (c) shows S as a function 
of time. As is expected with AFP-like schemes, population accumulates steadily in the 
desired state to extremely high fidelity. In keeping with the previous results, triplet readout 
goes to 1, and states in an equal superposition of singlet and triple asymptotically reach 0.5. 
The singlet state does not transfer population to the readout state, evinced by the observed 
S=0 across the singlet readout. The fact that there are still 'wiggles' in S suggest that some 
optimisation of the AFP scheme is) still desirable. 



20 After analysing the Digital Kane (DK) scheme [49], it has become clear that the spin 
readout scheme may be applicable to this scheme. The reasoning is that the triplet dot 
readout requires a three spatial eigenstates, which could be realised using donors to create 
the confining potential, or wells (for example in quantum dot structures). Our scheme 
could be comprised of a hybrid architecture, namely of donors and wells, for example 

25 where the Ir system is a two-donor system, and p is a quantum dot. Such an 
implementation would naturally arise in DK where qubits are defined by nuclear and 
electronic spins, and there are shuttling sites which are effective quantum dots at the 
Si/Si02 interface. It may therefore be possible to use this technique as a readout scheme for 
this alternative computing architecture. 

30 
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As stated before, the ideas presented in this work are highly general and should be 
applicable to any system where energy resolved states exist. We have concentrated our 
discussions on systems where the energy separation arises due to symmetric and anti- 
symmetric combinations of charge-position states, or to Cooper-Pair occupancies. 

It will be appreciated by persons skilled in the art that numerous variations and/or 
modifications may be made to the invention as shown in the specific embodiments without 
departing from the spirit or scope of the invention as broadly described. The present 





embodiments 


are, therefore, to be considered in all respects as illustrative and not 


10 


restrictive. 
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